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Abstract 

The space-like penguin contributions to branching ratios and CP asymmetries in 
charmless decays of B to two pseudoscalar mesons are studied using the next-to-leading 
order low energy effective Hamiltonian. Both the gluonic penguin and the electroweak 
penguin diagrams are considered. We find that the effects are significant. 
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1. Introduction 



Penguin diagrams play an important role in charmless B-decays and direct CP vio- 
lation [1,2]. But only time-like penguin diagrams were considered in the literature be- 
cause they can provide the necessary different strong phases for CP violation by different 
loop effects of the internal u and c quarks [1]. The contribution of space-like penguin 
diagrams is usually neglected assuming form factor suppression. This assumption for 
neglecting space-like penguin effects is used not only for gluonic penguins but also for 
electro- weak penguins [3]. But it does not lie on a solid ground because the space-like 
penguin amplitudes can be remarkably enhanced by the hadronic matrix elements involv- 
ing (V-A)(V+A) or (S+P)(S-P) currents [4]. Although space-like penguin diagrams can 
only provide an overall CP conserving phase due to final state interaction, it affects CP 
asymmetry by modifying the dispersive or absorptive parts of time-like penguin ampli- 
tudes. Obviously, it affects branching fractions too. In our recent paper [5], we illustrated 
the space-like penguin effects in CP asymmetries for the exclusive B-decays B~ — > K°ir~ 
and K°K~ using leading order Hamiltonian. In contrast to the naive expectation, the 
space-like penguin effects on CP asymmetries are found to be significant. In this letter 
we study space-like penguin effects in B to two pseudoscalar decays systematically. We 
concentrate on the charmless B decays because penguins play an important role in these 
decays. We use the next to leading order low energy effective Hamiltonian in order to 
consider both gluonic and electro-weak penguins. We can see later that the contribution 
of the electro-weak penguin is not negligible. The article is organized as following: In sec- 
tion II, We present the next to leading order effective Hamiltonian and the computation 
method. Section III devotes to the numerical results and corresponding discussions. 
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2. Effective hamiltonian and factorization 



approximation 

Following ref. [6], the next-to- leading order low energy effective Hamiltonian describ- 
ing AB = — 1 , AC = AU = transitions is given at the renormalization scale p = 0(m&) 



as 



n eff (AB = -i) = ^§ 



10 



]T v q q\cm + qicm + QkCM 

q=u,c I fc=3 



where the Wilson coefficient functions Cj(/x) (i=l,- • -,10) are calculated in renormalization 
group improved perturbation theory and include leading and next-to-leading order QCD 
corrections and leading order corrections in a. The CKM factors v q are defined as 



V* d V q b forbad transitions 



(2) 



V* s V q b for b — > s transitions. 

Here, we make use of the Wolfenstein parametrization[7] in which the CKM matrix can 
be written in terms of four parameters A, A, p and 77 in the following form: 



(3) 



A \ 3 A(p - irj) 
» = -A 1-|A 2 \ 2 A 

\ 3 A(l-p-ir]) -X 2 A 1 

The preferred values of the CKM Parameters are A = 0.22, A = 0.8, i] = 0.34, p = 
—0.12, which are obtained from the fit to experimental data[8]. The operators Q™, Q%, 
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Qzi ■ ■ ■ i Qio are given as the following forms: 



o u — 
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Q 7 = 


l(qb) v _ A J2e q ,(q , q') v+A 
<?' 


Q 8 = 


U^h)v-Aj2 e q'Wa) 
q' 


Q9 = 


l(qb) V - A ^e q ,(q'q , ) V - A 


Qw — 


U^)v-Aj2 e q'(yW*) 



(4) 



where Q" and Q2 are the current- current operators, and the current-current operators Ql 
and Q2 are obtained from and Q2 through the substitution of u — > c. Q3, . . . ,Qe are 
the QCD penguin operators, whereas Q7, • • • , Qio are the electroweak penguin operators. 
The quark g = d or s is for b — > (i or s transitions, respectively; q' is running over the quark 
flavours being active at the scale \i = 0{rrib) (q' G {u, d, c, s, b}); e q r are the corresponding 
quark charges; the indices a, (5 are SU(3) C color indices; (V ± A) refer to 7^(1 ± 75). 
It should be noted that the Hamiltonian (1) can be viewed as the generalization of the 
leading logarithmic Hamiltonians presented in [9,10]. 

Beyond the leading logarithmic approximation, the Wilson coefficient functions Cj(/i) 
depend both on the form of the operator basis (4) and on the renormalization scheme. 
Here, we use the renormalization scheme independent Wilson coefficient functions [11]: 



CM 



1 + Q »M f T + 



47T 



47T 



C(/x), 



(5) 



where r s and f e are obtained from one-loop matching conditions. Now, taking the 
QCD and electroweak one-loop level matrix elements of the operators Qi (Qi—Qi, Q2, 
Q3, ■ ■ ■ , Qw) into account through 



< Q T (/i) >=< Q T > 



• 



em ^ t 



Air 



m e (//) 



(6) 
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which define matrices m s (/i) and ra e (/i). In Eq. (5) and (6), C(/x), C(/i) and Q are all 
column vectors, where the vector Q are given by the operator basis Q iy and < Q > 
denote the tree level matrix elements of these operators. Combine Eq.(5) and (6), we 
obtain 



< Q T (^) • CM > 



<Q T >o 



1 + ^ (rhM - f s f + ^ (m e (/i) - f e f 



47T 



47T 



= < Q r > -C'(/i) 
where C'(//) are defined as 



CM 



(7) 



^9 



/^r /nr/ 

Ol, U 2 — L/ 2 , 

C3 — -Ps/3, C 4 = C4 + P s , 

C^ — Pg/S, Cg = 6*6 + -P s , 
C*7 + -fej C' 8 = Cg, 



(8) 



where P se are given by 



P s = %C 2 (ii) f-GKw) 



P e = ^(3C 1 + C 2 (/x 
G(m,q,n) = —A^dxx{l — x)ln 



f - G(m q ,q,n) 
m 2 — x(l — x)q 2 

V 2 



(9) 



here q — u, c, for numerical calculation, we take m u = 0.005GeV, m c = 1.35GeV, and q 2 
denotes the momentum transfer squared of the virtual gluons, photons and Z° appearing 
in the QCD and electroweak penguin matrix elements. For the details of this calculation, 
the reader is referred to ref. [12,13]. 

The renormalization scheme independent Wilson coefficient functions Ci(fj) at the 
scale pi = 0(m b ) are obtained by first calculating the Wilson coefficients at /i — 0(m w ) 
and then using the renormalization group equation to evolve them down to 0{m b ). We 



use in our analysis, a s (m z ) = 0.118, a(m z ) = 1/128[14] and m t = 174GeV[15] and the 
numerical values of the renormalization scheme independent Wilson Coefficients Cj(/i) at 
\i = 0{mb) are [13] 

ci = -0.313, c 2 = 1.150, c 3 = 0.017, 
c 4 = -0.037, c 5 = 0.010, c 6 = -0.046, 

(10) 

c 7 = -0.001 -a em , c 8 = 0.049 • a em , c 9 = -1.321 • a em , 
cio = 0.267 ■ a em . 

With the help of Eq. (7), the two-body decay amplitude < PP'\H eff (AB = -1)\B > 
can be expressed as linear combinations of < PP\Qi\B > . The hadronic matrix elements 
< PP \Qi\B > are evaluated using the factorization approximation [16]. It should be 
noted that this approach has already been used in the literature to analyze the QCD or 
electroweak time-like penguin contributions [12]. However, we go further in this letter by 
including the space-like penguin diagrams. As in [5,12], we also neglect W-annihilation or 
W-exchange diagram contributions in our present analysis which are commonly assumed 
to be form factor suppressed. 

Using the vacuum-saturation approximation, the decay amplitude < PP'\H e ff\B > 
can be factorized into a product of two current matrix elements < P\J^\0 > and < 
P'\J'\B > for the tree and time-like penguin diagrams (Fig.l), or the product of < 
pp'\ J M |0 > and < 0| J' \ B > for the space-like penguin diagrams (Fig. 2). In this work, the 
hadronic matrix elements are calculated in BSW method[16]. We define 

M ™ = < P\(QiQ2)v-a\0 >< P'\(bq 3 ) V -A\B >, 
< P\(qiq2)v-A\0 >< P'\{qzb)v- A \B >, 



11) 



or 



and 



S&w = < PP'\(QiQ2)v-a\0 >< 0\(bq 3 ) V - A \B >, 
= < PP'\(q iq2 ) v _ A \0 >< 0\(q 3 b) v _ A \B >, 
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(12) 



where denotes the hadronic matrix element in tree and time-like penguin diagram 

case, while S^f q2q3 denotes space-like penguin case. When the (V-A)(V+A) current are 
transformed into (S+P)(S-P) and further into (V-A)(V-A) ones using equation of motion 

for the time-like and space-like penguin amplitudes, there appear the terms which are 

i 1m\ 2m| 

proportional to — and — ; -, respectively. It q = q 

(m q + m q >)(m b - m q >) (m q - m q ,){m b + m q >) 

as in the decay modes: 

— > 7T _ 7r + , 7T 7r°, -r°t], n°if' , r/r], ifif , i]i], K° K° , ; (for b — > d transitions) 
— > K~K + , K°K°, ir°r], tc°7]' , r/rj, ifif , if' if'. (for b — > s transitions), 

2m 2 

the denominator of the factor — is zero. So, we can not use equation 

(m q - m ql ){m h + m q >) 

of motion to compute the amplitudes of these decays. We have to compute the matrix 
elements of (S+P)(S-P) operators directly. We shall discuss it elsewhere. 

As an example of how to factorize the decay amplitudes into the product of hadronic 
matrix elements, we give the result of < ir~ir°\H e ff\B~ > in the following, 

< 7i-7i°\H eff \B- > 



q=u,c 



efflu 

2M 2 

+ ° 3 - {m d + m u )(m u -m b ) ^ + ^ + ° 9 > M ^ + 



(13) 



M 2 3 3 1 o 

+{a 2 5« g - a 3 H ^ -(a 5 - a 7 /2) - -a 8 + -a w + -a 9 }M^J + 

m d [m d -m b ) 2 2 2 

+ {<■» + (m „ + ,^_ m „) ^ + «r) + «•}(«£.■" + «:r> 

where the term (S^~^° + ) is the contribution obtained from two space-like penguin 
diagrams, and the quark masses are taken as m d = O.OlGeV, m u = 0.005GeV, m s = 
0.175GeV, m b = 5.0GeV. a k is defined as 

a' 

aii-i = — r ^2ii 

a 2t = C^ + ^f, (i = l,2,3,4,5) 



The general expression for the one-body pseudoscalar matrix element of the axial-vector 



is 



<0\V^- A^\P(q) >=if P q li , 



(14) 



where q represents the momentum of the pseudoscalar meson, and fp is the decay con- 
stant. The two-body pseudoscalar-pseudoscalar matrix element of the vector current is[4, 
17] 

< P2(Q2)K - A.IPM >= U(q-)q + » + f-(q 2 -)q-», (is) 



where q± = qi ± q 2 , and the form factor f± is given by the monopole parametrization 

, / + (0) 



f-(d) 



1 — q 2 /m 2 



pole 



(16a) 
(16b) 



With equation (14~16), we obtain 



- if uu f B u «- {M 2_ 



M B - + 



(M B - + M n0 f - Ml 
(M B - + M n -) 2 - M 2 



(17) 



and 



S duu - 71 fB-fli M \ 



M n - -M n o 
B -' M w - +M w o 



D duu - J2JB-J+\ 1V1 B -) yr 



M.0 - Mrr 



(M w - + M n o) 2 - Mi- 



(M w - + M n0 ) 2 - Ml 



(18) 



M n - + M^o 

where the factors — come from the constituent of n° = -^(uu — dd). 

In order to give numerical results, we need to know the form factors. For the decay 
form factors like f+ n (M 2 ), etc., we can use BSW [16] method to calculate them. For 
the annihilation form factor /"(Q 2 ), we do not have reliable method to compute it. But 
at Q 2 = M B one is far from the Kti, 7T7t, r\K resonance region. So, for the charmless 
B decays, because the large energy release, we can use the form factor in its asymptotic 
form. For charmless B to two pseudoscalars decays, the asymptotic form factor predicted 
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by QCD [18] should be a resonable approximation. So we take f+(Q 2 ) = HQira s fg/Q 2 . 
Now we are in a position to give the numerical results. 



3. Results and discussions 



The decay width of a B-meson at rest decaying into two pseudoscalars is 

T(B - PP') = ±\< PP'\H eff \B > | 2 ^, (19) 

where 

= [(Ml - (Mp + M P 'f){Ml - (M P - M P 'f)]h 
m 2Mb 

is the momentum of the pseudoscalar meson P or P'. The corresponding branching ratios 
are given by 

Br(B -> PP') = T ( B ^ PP '\ (21) 

In our numerical calculation, we take[14] r£j = 4.27 x l(T 13 GeV, rf} t = 4.39 x l(T 13 GeV, 
and Tfj = 4.91 x l(T 13 GeV. 

In order to obtain the CP-violating parameter, the B-meson decay amplitude can be 
generally expressed as 

< PP'\H eff \B > = ^ Yl < C 'i <Qi> + C 2 <Ql> 

V ^ q=u,c 
10 

+ ]T C' k < Q k >) (22) 

fc=3 

= L V^- 

V ^ g=u,c 

With the help of eq.(22), one can get the CP- violating asymmetry parameter 

T(B -> PP') - T(P -> PP') 

(23) 



^ r(p -> pp') + r(P -> PP') 



2Jm( V ;)Jm(|) 



^ + ^(f) 2 + 2Pe(f)' 
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Since the branching ratios and CP asymmetries depend crucially on the parameter 
q 2 describing the momentum squared of the exchanged virtual particles appearing in the 
penguin matrix elements of Fig. 1 and 2, we should consider it in detail. Here, we use the 
same simple picture for two-body decays illustrated in Fig. lc and 2 as the one in ref [5]. 
With the simple physical picture presented in Ref. [5] , the average value of the momentum 
squared < q 2 > of the exchanged virtual particles can be given by 

< q 2 >= m 2 b + m 2 - 2m b E q , (24) 



where E q is determined from 



E q + JE* -m 2 q + m 2 , + - 4m 2 q + m 2 , = m b ; (25a) 



for the time-like penguin channels; or from 



E q + ^E 2 - m 2 + m 2 , = m b + m q > ; (25b) 

for the space-like penguin channels. When we factorize < Qk >o of hairpin diagrams 
illustrated in Fig. Id), we find < Q 3 > = — < Q5 >o, < Qa >o= — < Qs >q, and 
< Q 7 > = — < Qg >o- and hence the factor in Eq. 6: 



8tt 



jj- < Qs >o + < Qa >o < Q5 >o + < Qe >o 



Oi 

+— [< Q 7 > + < Q 9 > ] 

D7T 



(26) 



vanishes because of the cancelation. So, we do not need to consider the hairpin diagrams. 

The numerical results of the space-like penguin contributions to the branching ratios 
and CP- violating asymmetries are given in table 1. In the meantime, we also calculate the 
branching ratios and CP-violating asymmetries with only the tree and time-like penguin 
contributions for comparison. We also present the results with only tree and gluonic 
penguin contributions. All the parameters such as meson decay constants and form factors 
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needed in our calculation are taken as f n ± = 0.13GeV, f K = 0.160Gey[14], fffi = —fffl = 
U±/V2. f™ = /f = 0.092, /** = -0.105, /T = ff = 0.049, f°? = 0.12[4], f D = 0.23, 
f Ds = 0.281[19], f B = 1.5 x / ff± [20], f Bs = 0.206[21],and f^ n ~(0) = 0.29, f^ K ~(0) = 
0.32 [22], /+"^ u (0) = 0.307, (0) = 0.254, /f D (0) = 0.690 [16], /f r ""(0) = 0.335, 

/+•""'• (0) = 0.282 [23], /f aDs (0) = 0.648 [24]. 

From Table 1 we can see the following features: 

(i) For most of the charmless decays, penguin contributions are important. 

(ii) For B~ -> tt't], K'tt , K~r], K~rf, -> K°7r°, 7?% #V and B° s -> Tr ^ , the 
contribution of the electro-weak penguins are not negligible. 

(hi) The space-like penguin effects in B~ — > 7r~7r° are amazingly large. The correction 
to the branching ratio is more than 100%, while to the CP asymmetry is about an order of 
magnitude, actually, A cp ~ 0.77% with only time-like penguin, but A cp ~ —70.4% when 
including the space-like penguin. 

For B~ -> K-7i°, K-r], K~rf, KH~ , B° d -> ^-tt+, K%K°rf, and 5° -> 7r-fT+, 
Tr ^ , the space-like penguin contributions are also dominant. 

In B~ — > 7r~?7, 7T — 77', the space-like penguin contribution is zero. The reason is that 
there are two space-like penguin diagrams in each channel and the contributions of the 
two diagrams exactly cancel each other. 

In general, we can conclude that the space-like penguin effects are not negligible in 
most of the charmless two-pseudoscalar decays of the B mesons. The space-like penguins 
can affect not only CP asymmetries, but also decay branching ratios. Further investiga- 
tions are definitely needed. 
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Table 1. The branching ratios and the CP asymmetries, where the "Only Tree" means 
the branching ratios with only tree diagram contribution, "T-like" denotes the time-like 
penguin contributions, the "S-like" denotes the space-like penguin contributions, "QCD" 
means QCD penguin and tree diagrams contributions, and "QCD+EW" denotes full tree, 
QCD and EW (electro- weak) penguin contributions. 
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decay mode 


Br 




Only Tree 


Tree+T-like 


Tree+T-like+S-like 


Tree+T-like 


QCD 


QCD+EW 


QCD 


QCD+EW 


QCD 


QCD+E 1 




2.96 x 10~ 6 


2.90 x 1(T 6 


2.74 x 10~ 6 


1.30 x 10~ 6 


1.12 x lO" 6 


0.69% 


0.77% 


B u -»• 


2.19 x 10~ 6 


3.68 x 10- 6 


2.81 x 10~ 6 


3.68 x 10~ 6 


2.81 x lO" 6 


35.1% 


33.5% 


B- -> TT-l/ 


6.95 x 1(T 7 


5.84 x 1(T 6 


5.8 x 10~ 6 


5.84 x 10~ 6 


5.8 x 10~ 6 


16.0% 


16.1% 







4.16 x 10~ 7 


4.08 x 10~ 7 


3.79 x 10- 7 


3.72 x 10~ 7 


2.47% 


2.49% 




2.12 x 10~ 7 


2.93 x 1(T 6 


4.22 x 10~ 6 


5.29 x 10~ 6 


6.55 x 10~ 6 


-8.58% 


-6.15% 


B- -> 


1.57 x 1(T 7 


1.53 x 1(T 7 


1.84 x 10~ 7 


4.00 x 10~ 8 


7.13 x 10~ 8 


6.17% 


4.38% 




4.98 x 1(T 8 


8.19 x lO" 6 


7.74 x 10~ 6 


8.03 x 10- 6 


7.59 x 10- 6 


-3.35% 


-3.53% 







4.95 x 1(T 6 


4.86 x 10~ 6 


7.10 x 10~ 6 


6.97 x 10~ 6 


-0.18% 


-0.18% 




3.77 x 10- 7 


5.78 x 10- 6 


5.96 x 10~ 6 


1.04 x 10~ 5 


1.06 x lO" 5 


-8.23% 


-8.0% 




6.76 x 1(T 10 


2.46 x 1(T 6 


1.56 x 10~ 6 


3.64 x 10~ 6 


2.74 x lO" 6 


0.44% 


0.73% 


B° d -> ^ 


6.72 x 1(T 10 


4.25 x 1(T 8 


1.97 x 10~ 9 


6.42 x 10~ 8 


2.33 x 10~ 8 


1.40% 


47.0% 


£° - K°r/ 


1.87 x 10~ 10 


7.97 x 10- 6 


7.42 x 10~ 6 


7.64 x 10~ 6 


7.09 x 10" 6 


-0.42% 


-0.44% 




3.51 x 1(T 6 
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Figure Captions 



Fig. 1. Quark diagrams for a B meson decaying into two light pseudoscalar mesons P 
and P' through the tree process b — > u{uq): a) the internal ^-emission diagram, (b) the 
external jy-emission diagram; and the time-like penguin diagram process b — > q(q'q')' c) 
the time-like pure penguin diagram, and d) the time-like hairpin diagram. The subscripts 
"s" denote "spectator" . The dark dot stands for the contraction of the W-loop. 
Fig. 2. Quark diagrams for a B meson decaying into two light pseudoscalar mesons P 
and P' through the space-like penguin process (bq r ) — > (qq'). The subscripts "v" denote 
"vacuum". The dark dot stands for the contraction of the W-loop. 
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Fig. 2. 
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